Evidence of the interplay of genetics and culture in Ethiopia.

The rich linguistic, ethnic and cultural diversity of Ethiopia provides an unprecedented opportunity to understand the level to which cultural factors correlate with-and shape-genetic structure in human populations. Using primarily new genetic variation data covering 1,214 Ethiopians representing 68 different ethnic groups, together with information on individuals' birthplaces, linguistic/religious practices and 31 cultural practices, we disentangle the effects of geographic distance, elevation, and social factors on the genetic structure of Ethiopians today. We provide evidence of associations between social behaviours and genetic differences among present-day peoples. We show that genetic similarity is broadly associated with linguistic affiliation, but also identify pronounced genetic similarity among groups from disparate language classifications that may in part be attributable to recent intermixing. We also illustrate how groups reporting the same culture traits are more genetically similar on average and show evidence of recent intermixing, suggesting that shared cultural traits may promote admixture. In addition to providing insights into the genetic structure and history of Ethiopia, we identify the most important cultural and geographic predictors of genetic differentiation and provide a resource for designing sampling protocols for future genetic studies involving Ethiopians.

Contrasting exome constancy and regulatory region variation in the gene encoding CYP3A4: an examination of the extent and potential implications.

OBJECTIVE: CYP3A4 expression varies up to 100-fold among individuals, and, to date, genetic causes remain elusive. As a major drug-metabolizing enzyme, elucidation of such genetic causes would increase the potential for introducing personalized dose adjustment of therapies involving CYP3A4 drug substrates. The foetal CYP3A isoform, CYP3A7, is reported to be expressed in ∼10% of European adults and may thus contribute towards the metabolism of endogenous substances and CYP3A drug substrates. However, little is known about the distribution of the variant expressed in the adult. METHODS: We resequenced the exons, flanking introns, regulatory elements and 3'UTR of CYP3A4 in five Ethiopian populations and incorporated data from the 1000 Genomes Project. Using bioinformatic analysis, we assessed likely consequences of observed CYP3A4 genomic variation. We also conducted the first extensive geographic survey of alleles associated with adult expression of CYP3A7 - that is, CYP3A7*1B and CYP3A7*1C. RESULTS AND CONCLUSION: Ethiopia contained 60 CYP3A4 variants (26 novel) and more variants (>1%) than all non-African populations combined. No nonsynonymous mutation was found in the homozygous form or at more than 2.8% in any population. Seventy-nine per cent of haplotypes contained 3'UTR and/or regulatory region variation with striking pairwise population differentiation, highlighting the potential for interethnic variation in CYP3A4 expression. Conversely, coding region variation showed that significant interethnic variation is unlikely at the protein level. CYP3A7*1C was found at up to 17.5% in North African populations and in significant linkage disequilibrium with CYP3A5*3, indicating that adult expression of the foetal isoform is likely to be accompanied by reduced or null expression of CYP3A5.

Evidence for a Common Origin of Blacksmiths and Cultivators in the Ethiopian Ari within the Last 4500 Years: Lessons for Clustering-Based Inference.

The Ari peoples of Ethiopia are comprised of different occupational groups that can be distinguished genetically, with Ari Cultivators and the socially marginalised Ari Blacksmiths recently shown to have a similar level of genetic differentiation between them (FST ≈ 0.023 - 0.04) as that observed among multiple ethnic groups sampled throughout Ethiopia. Anthropologists have proposed two competing theories to explain the origins of the Ari Blacksmiths as (i) remnants of a population that inhabited Ethiopia prior to the arrival of agriculturists (e.g. Cultivators), or (ii) relatively recently related to the Cultivators but presently marginalized in the community due to their trade. Two recent studies by different groups analysed genome-wide DNA from samples of Ari Blacksmiths and Cultivators and suggested that genetic patterns between the two groups were more consistent with model (i) and subsequent assimilation of the indigenous peoples into the expanding agriculturalist community. We analysed the same samples using approaches designed to attenuate signals of genetic differentiation that are attributable to allelic drift within a population. By doing so, we provide evidence that the genetic differences between Ari Blacksmiths and Cultivators can be entirely explained by bottleneck effects consistent with hypothesis (ii). This finding serves as both a cautionary tale about interpreting results from unsupervised clustering algorithms, and suggests that social constructions are contributing directly to genetic differentiation over a relatively short time period among previously genetically similar groups.

Genetic signatures reveal high-altitude adaptation in a set of ethiopian populations.

The Tibetan and Andean Plateaus and Ethiopian highlands are the largest regions to have long-term high-altitude residents. Such populations are exposed to lower barometric pressures and hence atmospheric partial pressures of oxygen. Such "hypobaric hypoxia" may limit physical functional capacity, reproductive health, and even survival. As such, selection of genetic variants advantageous to hypoxic adaptation is likely to have occurred. Identifying signatures of such selection is likely to help understanding of hypoxic adaptive processes. Here, we seek evidence of such positive selection using five Ethiopian populations, three of which are from high-altitude areas in Ethiopia. As these populations may have been recipients of Eurasian gene flow, we correct for this admixture. Using single-nucleotide polymorphism genotype data from multiple populations, we find the strongest signal of selection in BHLHE41 (also known as DEC2 or SHARP1). Remarkably, a major role of this gene is regulation of the same hypoxia response pathway on which selection has most strikingly been observed in both Tibetan and Andean populations. Because it is also an important player in the circadian rhythm pathway, BHLHE41 might also provide insights into the mechanisms underlying the recognized impacts of hypoxia on the circadian clock. These results support the view that Ethiopian, Andean, and Tibetan populations living at high altitude have adapted to hypoxia differently, with convergent evolution affecting different genes from the same pathway.

Molecular diversity and population structure at the Cytochrome P450 3A5 gene in Africa.

BACKGROUND: Cytochrome P450 3A5 (CYP3A5) is an enzyme involved in the metabolism of many therapeutic drugs. CYP3A5 expression levels vary between individuals and populations, and this contributes to adverse clinical outcomes. Variable expression is largely attributed to four alleles, CYP3A5*1 (expresser allele); CYP3A5*3 (rs776746), CYP3A5*6 (rs10264272) and CYP3A5*7 (rs41303343) (low/non-expresser alleles). Little is known about CYP3A5 variability in Africa, a region with considerable genetic diversity. Here we used a multi-disciplinary approach to characterize CYP3A5 variation in geographically and ethnically diverse populations from in and around Africa, and infer the evolutionary processes that have shaped patterns of diversity in this gene. We genotyped 2538 individuals from 36 diverse populations in and around Africa for common low/non-expresser CYP3A5 alleles, and re-sequenced the CYP3A5 gene in five Ethiopian ethnic groups. We estimated the ages of low/non-expresser CYP3A5 alleles using a linked microsatellite and assuming a step-wise mutation model of evolution. Finally, we examined a hypothesis that CYP3A5 is important in salt retention adaptation by performing correlations with ecological data relating to aridity for the present day, 10,000 and 50,000 years ago. RESULTS: We estimate that ~43% of individuals within our African dataset express CYP3A5, which is lower than previous independent estimates for the region. We found significant intra-African variability in CYP3A5 expression phenotypes. Within Africa the highest frequencies of high-activity alleles were observed in equatorial and Niger-Congo speaking populations. Ethiopian allele frequencies were intermediate between those of other sub-Saharan African and non-African groups. Re-sequencing of CYP3A5 identified few additional variants likely to affect CYP3A5 expression. We estimate the ages of CYP3A5*3 as ~76,400 years and CYP3A5*6 as ~218,400 years. Finally we report that global CYP3A5 expression levels correlated significantly with aridity measures for 10,000 [Spearmann's Rho= -0.465, p=0.004] and 50,000 years ago [Spearmann's Rho= -0.379, p=0.02]. CONCLUSIONS: Significant intra-African diversity at the CYP3A5 gene is likely to contribute to multiple pharmacogenetic profiles across the continent. Significant correlations between CYP3A5 expression phenotypes and aridity data are consistent with a hypothesis that the enzyme is important in salt-retention adaptation.

The frequency of an IL-18-associated haplotype in Africans.

Variation within the gene for the proinflammatory cytokine interleukin (IL)-18 has been associated with inter-individual differences in levels of free protein and disease risk. We investigated the frequency of function-associated IL18 gene haplotypes in an extensive sample (n=2357) of African populations from across the continent. A previously identified five tagging SNP (single-nucleotide polymorphism) haplotype (here designated hGTATA), known to be associated with lower levels of IL-18, was observed at a frequency of 27% in a British population of recent European ancestry, but was found at low frequency (<8%) or completely absent in African populations. Potentially protective variants may, as a consequence, be found at low frequency in African individuals and may confer a difference in disease risk.

Ethiopian genetic diversity reveals linguistic stratification and complex influences on the Ethiopian gene pool.

Humans and their ancestors have traversed the Ethiopian landscape for millions of years, and present-day Ethiopians show great cultural, linguistic, and historical diversity, which makes them essential for understanding African variability and human origins. We genotyped 235 individuals from ten Ethiopian and two neighboring (South Sudanese and Somali) populations on an Illumina Omni 1M chip. Genotypes were compared with published data from several African and non-African populations. Principal-component and STRUCTURE-like analyses confirmed substantial genetic diversity both within and between populations, and revealed a match between genetic data and linguistic affiliation. Using comparisons with African and non-African reference samples in 40-SNP genomic windows, we identified "African" and "non-African" haplotypic components for each Ethiopian individual. The non-African component, which includes the SLC24A5 allele associated with light skin pigmentation in Europeans, may represent gene flow into Africa, which we estimate to have occurred ~3 thousand years ago (kya). The non-African component was found to be more similar to populations inhabiting the Levant rather than the Arabian Peninsula, but the principal route for the expansion out of Africa ~60 kya remains unresolved. Linkage-disequilibrium decay with genomic distance was less rapid in both the whole genome and the African component than in southern African samples, suggesting a less ancient history for Ethiopian populations.

Prevalence of clinically relevant UGT1A alleles and haplotypes in African populations.

Variation of a short (TA)(n) repeat sequence (rs8175347) covering the TATA box of UGT1A1 (UDP-glucuronosyltransferase1A1) is associated with hyperbilirubinaemia (Gilbert's syndrome) and adverse drug reactions, and is used for dosage advice for irinotecan. Several reports indicate that the low-activity (risk) alleles ((TA)(7) and (TA)(8) )) are very frequent in Africans but the patterns of association with other variants in the UGT1A gene complex that may modulate these responses are not well known. rs8175347 and two other clinically relevant UGT1A variants (rs11692021 and rs10929302) were assayed in 2616 people from Europe and Africa. Low-activity (TA)(n) alleles frequencies were highest in equatorial Africa, (TA)(7,) being the most common in Cameroon, Ghana, southern Sudan, and in Ethiopian Anuak. Haplotypic diversity was also greatest in equatorial Africa, but in Ethiopia was very variable across ethnic groups. Resequencing of the promoter of a sample subset revealed no novel variations, but rs34547608 and rs887829 were typed and shown to be tightly associated with (TA)(n) . Our results illustrate the need for investigation of the effect of UGT1A variants other than (TA)(n) on the risk of irinotecan toxicity, as well as hyperbilirubinaemia due to hemolytic anaemia or human immunodeficiency virus protease inhibitors, so that appropriate pharmacogenetic advice can be given.

CYP1A2 is more variable than previously thought: a genomic biography of the gene behind the human drug-metabolizing enzyme.

BACKGROUND AND OBJECTIVES: CYP1A2 metabolizes various drugs, endogenous compounds and procarcinogens. As human genetic diversity has been reported to decrease with distance from Ethiopia, we resequenced CYP1A2 in five Ethiopian ethnic groups representing a rough northeast to southwest transect across Ethiopia to establish: (i) what variation exists in comparison with what is already known globally and (ii) what CYP1A2 pharmacogenetic profiles may be present as several CYP1A2-metabolized drugs are administered to Ethiopians. RESULTS AND CONCLUSIONS: We found 49 different variable sites (30 of which are novel), nine nonsynonymous changes (seven of which are novel), one synonymous change and 55 different haplotypes, only three of which are previously reported. When haplotypes were constructed using only nonsynonymous polymorphisms to restrict haplotypes to those most likely to affect enzyme structure/function, 10 haplotypes were identified (seven contain previously unidentified nonsynonymous variants and four are predicted to alter the enzyme structure/function). Most individuals have at least one copy of the ancestral haplotype. Comparing these data with those from publically available databases, Ethiopian groups display twice the variation seen in all other populations combined (gene diversity using nonsynonymous variants): Ethiopia=0.17±0.02, other populations=0.08±0.03. Across the entire gene, Ethiopia also evidences all common variation found on a global scale. We provide evidence of weak purifying selection acting on CYP1A2 and show that the time to most recent common ancestor, calculated using variation in a nearby microsatellite, places several variants into a period predating the expansion of modern humans out of Africa less than 100,000 years ago.

Missense mutations in the APOL1 gene are highly associated with end stage kidney disease risk previously attributed to the MYH9 gene.

MYH9 has been proposed as a major genetic risk locus for a spectrum of nondiabetic end stage kidney disease (ESKD). We use recently released sequences from the 1000 Genomes Project to identify two western African-specific missense mutations (S342G and I384M) in the neighboring APOL1 gene, and demonstrate that these are more strongly associated with ESKD than previously reported MYH9 variants. The APOL1 gene product, apolipoprotein L-1, has been studied for its roles in trypanosomal lysis, autophagic cell death, lipid metabolism, as well as vascular and other biological activities. We also show that the distribution of these newly identified APOL1 risk variants in African populations is consistent with the pattern of African ancestry ESKD risk previously attributed to MYH9.Mapping by admixture linkage disequilibrium (MALD) localized an interval on chromosome 22, in a region that includes the MYH9 gene, which was shown to contain African ancestry risk variants associated with certain forms of ESKD (Kao et al. 2008; Kopp et al. 2008). MYH9 encodes nonmuscle myosin heavy chain IIa, a major cytoskeletal nanomotor protein expressed in many cell types, including podocyte cells of the renal glomerulus. Moreover, 39 different coding region mutations in MYH9 have been identified in patients with a group of rare syndromes, collectively termed the Giant Platelet Syndromes, with clear autosomal dominant inheritance, and various clinical manifestations, sometimes also including glomerular pathology and chronic kidney disease (Kopp 2010; Sekine et al. 2010). Accordingly, MYH9 was further explored in these studies as the leading candidate gene responsible for the MALD signal. Dense mapping of MYH9 identified individual single nucleotide polymorphisms (SNPs) and sets of such SNPs grouped as haplotypes that were found to be highly associated with a large and important group of ESKD risk phenotypes, which as a consequence were designated as MYH9-associated nephropathies (Bostrom and Freedman 2010). These included HIV-associated nephropathy (HIVAN), primary nonmonogenic forms of focal segmental glomerulosclerosis, and hypertension affiliated chronic kidney disease not attributed to other etiologies (Bostrom and Freedman 2010). The MYH9 SNP and haplotype associations observed with these forms of ESKD yielded the largest odds ratios (OR) reported to date for the association of common variants with common disease risk (Winkler et al. 2010). Two specific MYH9 variants (rs5750250 of S-haplotype and rs11912763 of F-haplotype) were designated as most strongly predictive on the basis of Receiver Operating Characteristic analysis (Nelson et al. 2010). These MYH9 association studies were then also extended to earlier stage and related kidney disease phenotypes and to population groups with varying degrees of recent African ancestry admixture (Behar et al. 2010; Freedman et al. 2009a, b; Nelson et al. 2010), and led to the expectation of finding a functional African ancestry causative variant within MYH9. However, despite intensive efforts including re-sequencing of the MYH9 gene no suggested functional mutation has been identified (Nelson et al. 2010; Winkler et al. 2010). This led us to re-examine the interval surrounding MYH9 and to the detection of novel missense mutations with predicted functional effects in the neighboring APOL1 gene, which are significantly more associated with ESKD than all previously reported SNPs in MYH9.

Multiple rare variants as a cause of a common phenotype: several different lactase persistence associated alleles in a single ethnic group.

Persistence of intestinal lactase into adulthood allows humans to use milk from other mammals as a source of food and water. This genetic trait has arisen by convergent evolution and the derived alleles of at least three different single nucleotide polymorphisms (-13910C>T, -13915T>G, -14010G>C) are associated with lactase persistence in different populations. Each allele occurs on an extended haplotype, consistent with positive directional selection. The SNPs are located in an 'enhancer' sequence in an intron of a neighboring gene (MCM6) and modulate lactase transcription in vitro. However, a number of lactase persistent individuals carry none of these alleles, but other low-frequency single nucleotide polymorphisms have been observed in the same region. Here we examine a cohort of 107 milk-drinking Somali camel-herders from Ethiopia. Eight polymorphic sites are identified in the enhancer. -13915*G and -13907*G (a previously reported candidate) are each significantly associated with lactase persistence. A new allele, -14009*G, has borderline association with lactase persistence, but loses significance after correction for multiple testing. Sequence diversity of the enhancer is significantly higher in the lactase persistent members of this and a second cohort compared with non-persistent members of the two groups (P = 7.7 x 10(-9) and 1.0 x 10(-3)). By comparing other loci, we show that this difference is not due to population sub-structure, demonstrating that increased diversity can accompany selection. This contrasts with the well-documented observation that positive selection decreases diversity by driving up the frequency of a single advantageous allele, and has implications for association studies.

The potentially deleterious functional variant flavin-containing monooxygenase 2*1 is at high frequency throughout sub-Saharan Africa.

BACKGROUND: The drug-metabolizing enzyme flavin-containing monooxygenase 2 (FMO2) is the predominant FMO isoform present in the lung of most mammals, including non-human primates. All Europeans and Asians tested have been shown to be homozygous for a non-functional variant, FMO2*2A, which contains a premature stop codon due to a single-nucleotide change in exon 9 (g.23238C>T). The ancestral allele, FMO2*1, encodes a functionally active protein and has been found in African-Americans (26%) and Hispanics (2% to 7%). Possessing this variant increases the risk of pulmonary toxicity when exposed to thioureas, a widely used class of industrial compounds. FMO2 may also be involved in the metabolism of drugs that are used to treat diseases that are prevalent in Africa. RESULTS AND CONCLUSION: We conducted a survey of g.23238C>T variation across Africa that revealed that the distribution of this SNP is relatively homogeneous across sub-Saharan Africa, with approximately one third of individuals possessing at least one FMO2*1 allele, though in some populations the incidence of these individuals approached 50%. Thus many sub-Saharan Africans may be at substantially increased health risk when encountering thiourea-containing substrates of FMO2. Analysis of HapMap data with the Long-Range Haplotype test found no evidence for positive selection of either 23238C>T allele and maximum-likelihood coalescent analysis indicated that this mutation occurred some 500,000 years before present. This study demonstrates the value of performing genetic surveys in Africa, a continent in which human genetic diversity is thought to be greatest, but where studies of the distribution of this diversity are few.

A novel polymorphism associated with lactose tolerance in Africa: multiple causes for lactase persistence?

Persistence or non-persistence of lactase expression into adult life is a polymorphic trait that has been attributed to a single nucleotide polymorphism (C-13910T) in an enhancer element 13.9 kb upstream of the lactase gene (LCT). The -13910*T allele occurs at very high frequency in northern Europeans as part of a very long haplotype (known as A), and promotes binding of the transcription factor Oct-1. However, -13910*T is at very low frequency in many African milk drinking pastoralist groups where lactase persistence phenotype has been reported at high frequency. We report here for the first time, a cohort study of lactose digester and non-digester Sudanese volunteers and show there is no association of -13910*T or the A haplotype with lactase persistence. We support this finding with new genotype/phenotype frequency comparisons in pastoralist groups of eastern African and Middle Eastern origin. Resequencing revealed three new SNPs in close proximity to -13910*T, two of which are within the Oct-1 binding site. The most frequent of these (-13915*G) is associated with lactose tolerance in the cohort study, providing evidence for a cis-acting effect. Despite its location, -13915*G abolishes, rather than enhances Oct-1 binding, indicating that this particular interaction is unlikely to be involved in lactase persistence. This study reveals the complexity of this phenotypic polymorphism and highlights the limitations of C-13910T as a diagnostic test for lactase persistence status, at least for people with non-European ancestry.

The T allele of a single-nucleotide polymorphism 13.9 kb upstream of the lactase gene (LCT) (C-13.9kbT) does not predict or cause the lactase-persistence phenotype in Africans.

The ability to digest the milk sugar lactose as an adult (lactase persistence) is a variable genetic trait in human populations. The lactase-persistence phenotype is found at low frequencies in the majority of populations in sub-Saharan Africa that have been tested, but, in some populations, particularly pastoral groups, it is significantly more frequent. Recently, a CT polymorphism located 13.9 kb upstream of exon 1 of the lactase gene (LCT) was shown in a Finnish population to be closely associated with the lactase-persistence phenotype (Enattah et al. 2002). We typed this polymorphism in 1,671 individuals from 20 distinct cultural groups in seven African countries. It was possible to match seven of the groups tested with groups from the literature for whom phenotypic information is available. In five of these groups, the published frequencies of lactase persistence are >/=25%. We found the T allele to be so rare that it cannot explain the frequency of the lactase-persistence phenotype throughout Africa. By use of a statistical procedure to take phenotyping and sampling errors into account, the T-allele frequency was shown to be significantly different from that predicted in five of the African groups. Only the Fulbe and Hausa from Cameroon possessed the T allele at a level consistent with phenotypic observations (as well as an Irish sample used for comparison). We conclude that the C-13.9kbT polymorphism is not a predictor of lactase persistence in sub-Saharan Africans. We also present Y-chromosome data that are consistent with previously reported evidence for a back-migration event into Cameroon, and we comment on the implications for the introgression of the -13.9kb*T allele.

Founding mothers of Jewish communities: geographically separated Jewish groups were independently founded by very few female ancestors.

We have analyzed the maternally inherited mitochondrial DNA from each of nine geographically separated Jewish groups, eight non-Jewish host populations, and an Israeli Arab/Palestinian population, and we have compared the differences found in Jews and non-Jews with those found using Y-chromosome data that were obtained, in most cases, from the same population samples. The results suggest that most Jewish communities were founded by relatively few women, that the founding process was independent in different geographic areas, and that subsequent genetic input from surrounding populations was limited on the female side. In sharp contrast to this, the paternally inherited Y chromosome shows diversity similar to that of neighboring populations and shows no evidence of founder effects. These sex-specific differences demonstrate an important role for culture in shaping patterns of genetic variation and are likely to have significant epidemiological implications for studies involving these populations. We illustrate this by presenting data from a panel of X-chromosome microsatellites, which indicates that, in the case of the Georgian Jews, the female-specific founder event appears to have resulted in elevated levels of linkage disequilibrium.